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Abstract 

Fermat's principle applied to a flat metric in the plane yields the phase 
^ ' of a Bessel function in the periodic domain for a constant index of refrac- 

tion. Gravitational forces cause the index of refraction to vary and lead to 
a modified phase of the Bessel function. A distinction is made between the 
forces that cause acceleration: the gravitational force affects the optical 
properties of the medium whereas the centrifugal force does not, the latter 
• being built into the phase of oscillations of the Bessel function. The time 

delay in radar echoes from planets is determined from Fermat's principle 
^ where the velocity of light is the phase velocity and the index of refraction 

varies on account of the gravitational potential. The deflection of light by 
a massive body is shown to be produced by a quadrupole interaction, and 
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Or 



■ the perihelion shift requires both the gravitational potential, producing a 



closed orbit, and the quadrupole, causing the perihelion to rotate. 



1 Introduction 

X. 

General relativity is based on the notion that gravity, rather than being a force 
between masses, is the curvature of spacetime. The source of curvature is mass 
itself, like the source of an electric field is charge. Just as free particles follow 
straight lines in flat spacetime, free particles follow geodesies in the curved 
spacetime of a gravitational field. The calculations of the red shift, the time 
delay in radar echoes from planets, the bending of light, the perihelion shift, and 
the geodesic effect all support general relativity. 

Sometimes new insight can be gained in looking at old results from a new 
perspective. The old quantum theory of multiply-periodic central motion was 
superseded by wave mechanics. Old quantum theory attempted to apply celes- 
tial mechanics to the world of atoms with the addition of quantum restrictions to 
the motion Bo rn 1925] . In this article we will apply it to diffraction phenomena 
in the short-wavelength limit that involve the relativistic effects of time delay in 
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radar sounding, the deflection of light, and the advance of the perihelion — three 
of the so-called five tests of general relativity. Moreover, these relativistic effects 
are both static so that they can be investigated by studying only the spatial 
component of the metric. The trajectory of a light ray in a static gravitational 
field is determined in the same way as in an inhomogeneous refractive medium 
Fock 1966 . Format's principle of least time relates the length and orientation 
of a light ray to the time for light to propagate along the ray path. The analogy 
between the index of refraction and the square root of twice the difference be- 
tween the total and potential energies is also well known. Yet, this can provide 
a way of distinguishing between the centrifugal and gravitational fields which 
cause acceleration, and the form of the trajectory of the motion. 

If we take a flat space time metric in the plane, and consider a constant index 
of refraction, we will show that Fermat's principle yields precisely the phase of 
the oscillations of the Bessel function of the first kind in the periodic domain in 
the asymptotic, short- wavelength, limit. This is the same as the WKB result, 
and it allows us to associate a wave phenomenon with a geodesic trajectory. The 
only potential appearing is the repulsive centrifugal potential, and the trajectory 
is necessarily open. In contrast to general relativity, where gravity is not con- 
sidered as a force in the conventional sense, but, rather, built into the curvature 
of spacetime, the centrifugal force is built into the phase of the Bessel function 
in the periodic domain, where the trajectory consists of straight line segments 
and arc lengths on a caustic circle whose radius is determined by the magnitude 
of the angular momentum. Whereas all forces causing acceleration are on the 
same footing in general relativity, the gravitational force, and not centrifugal 
force, has to be introduced through a space varying index of refraction in an 
inhomogeneous refractive medium. The effect of gravity is to make the medium 
optically more dense in the vicinity of a massive body, while the centrifugal force 
has no effect on the optical properties of the medium. Centrifugal, Coriolis and 
gravitational forces are usually considered to be fictitious insofar as they can be 
transformed away by a change of frame. The centrifugal and Coriolis forces can 
be transformed away by changing to a nonrotating frame, while the fictitious 
force of gravity is transformed away by changing from a nonfreely falling to a 
freely falling frame. In Newtonian theory one insists on a nonrotating frame, 
but not a freely falling one. Here, we appreciate that the force of gravity affects 
the optical properties of the medium whereas the centrifugal force determines 
the radius of the caustic circle and the diffraction pattern in its neighborhood. 

The time delay in radar sounding, the bending of light in the vicinity of a 
massive body and the advance of the perihelion will be treated from this view- 
point. The introduction of the Newtonian potential will cause a modification 
of the phase of the Bessel function and yield both periodic and aperiodic orbits 
depending upon whether the total energy is negative or positive, respectively. 
The time delay in radar sounding is a double relativistic effect that is caused by 
the slowing down of clocks in the vicinity of a massive body, accounted for by 
the fact that light travels at the phase velocity, and the shrinking of measuring 
rods, accounted for by a line element that is magnified by the index of refraction. 
Since an index of refraction that takes into account the Newtonian gravitational 
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potential only gives one half the relativistic value it cannot be the cause, or at 
least the entire cause, of the bending of light. Rather, if we take the next higher 
interaction into account that causes mass to accelerate we do find the actual 
general relativistic value without the Newtonian gravitational potential. This 
implies that, like gravitational radiation, the interaction between a light ray and 
a massive body is predominantly quadrupole. 

In contrast, the advance of the perihelion requires both gravitational po- 
tential and the quadrupole interaction. In general relativity the quadrupole 
appears as a relativistic correction to the square of the transverse velocity in 
the conservation of energy. The gravitational potential is responsible for the 
closed elliptical orbit, while the quadrupole causes the perihelion to slowly ro- 
tate producing a rosette orbit. A dipole moment would have been sufficient to 
cause the advance of the perihelion |Born 1925] , but, since there is conservation 
of momentum, the center of mass of the system cannot accelerate, and so the 
mass dipole moment cannot accelerate. 



2 Metric and Phase 

Fermat's principle states that the ray path connecting two arbitrary points 
makes the optical path length 

I = J r]V2Tdt (1) 

stationary, where r\ is the index of refraction and T is the kinetic energy per 
unit mass. As a first application of Fermat's principle we calculate the time 
delay in radar sounding. In general relativity, this time delay is predicted by 
the Schwarzschild metric evaluated on the null geodesic when all angular depen- 
dencies are ignored. If a light signal is sent from the earth, located along the 
x-axis at —xe, to Venus, which is located behind the sun at xy, the light ray 
will be bent as it passes the gravitational field of the sun. Clocks will be slowed 
down and the time it takes the ray to bounce off Venus' surface and return to 
earth will be longer than if the sun were not present. 

The simplest mechanical analog of the index of refraction is r\ — — 2U/c 2 , 
where U is the potential energy per unit mass, and c is the speed of light in a 
vacuum. Since the gravitational field of the sun makes the medium optically 
more dense, U is identified as the gravitational potential —GM/r, where G is 
the Newtonian gravitational constant, M is the sun's mass, r = y/ R 2 + x 2 is 
the distance from the center of the sun to Venus, and R is the radius of the sun. 

Now, according to Fermat's principle, the propagation time t along a ray 
connecting the two endpoints —xe and xy is given by 

t = r vM d X = r VTTWr dx 

J-xb u ( r ) '-xe u(r) 

where 1Z := 2GM/c 2 is the Schwarzschild radius. The slowing down of clocks 
in a gravitational field will result in an apparent reduction in the speed of light. 
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Light will therefore travel at the phase velocity u(r) = c/rj(r), rather than c, as 
it would in vacuum. Consequently, the travel time will be 

K f Xv dx K f Xv dx 

t + At = t N + — — = t N + — 



c J- XE r c J_ XE Vi? 2 + x 2 

where tjsr — (xy + xe)/c is the Newtonian travel time. The second term is half 
the time delay for a signal to bounce off Venus and return to the earth. Fermat's 
principle thus predicts a time dilatation of 



a* = i — ^ m M (B) ln ( W| = 2 . 4 x 10 - 4 Bj 



c 



-xe 



J V c J \ r 2 J 



where the square roots have been expanded to lowest order using the fact that 
R <C xy, xe- 

This simple application of Fermat's principle gives the same result as general 
relativity, 1 and which is to within 3 % of experimental uncertainty |Shapiro~1 968 . 
In the following, we will use gravitational units where c = G = 1, and the more 
general expression for the index of refraction 



= V-A-2U, (2) 

where the dimensionless constant —A represents the negative of twice the total 
energy, allowing for both positive and negative values. 

In a rotating coordinate system in the Euclidean plane (9 = ir/2) 

2T = r 2 +r 2 <j) 2 . (3) 

Fermat's principle will determine the phase of the wave function tp that is a 
solution of Helmholtz's equation, 

(A + k 2 ) i) = 0, 

in the short-wavelength, or high-wave number k, limit, where A stands for the 
Laplacian. The solution to Helmholtz's equation can be written as 

^( r ) = ^ e ±^)(l + 0( K - 1 )), (4) 

where A is a constant term in the amplitude, ±(k5 i — •&) are phases of the 
incoming and outgoing waves, and $ is an integration constant that will be 

1 In a simplified demonstration of the time delay of radar sounding caused by a massive body 
JSexl an d Scxl 1979 , t = J dx/c e g, where c e g = (1 — TZ/2r) 2 c ss (1 — lZ/r)c, is an 'effective 
speed' that supposedly accounts for both time dilatation and the shrinking of measuring rods 
in a gravitational field. The final expression for t is valid to first order in TZ/r. Rather, in our 
presentation, the expression for t is exact, and it is a result of two factors: the slowing down of 
clocks, as expressed by the fact that light rays propagate at the phase velocity, c/rj(r), and the 
warping of measuring sticks, as accounted for by the stretching of the line element, r](r)dx. If 
the effective velocity c e g were to be identified as the phase velocity, then this would account 
for only half of the effect. 
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I = 



determined by matching conditions between the periodic and exponential do- 
mains [vid., l|15l) below]. The prime denotes differentiation with respect to r. 
The function S is referred to as the eikonal, and it is the object of our study. 

Introducing Q into Q implies that the ray path connecting two arbitrary 
points makes the optical length 

J 7]\Jr 2 + r 2 4> 2 dt = J r\\J\ + r 2 (j)' 2 dr (5) 

stationary. Observing that <fi is a cyclic coordinate, and calling A the integrand 
of JSJ , we immediately obtain a first integral of the motion 

<9A T]r 2 (/)' 

— — — — = = r a = constant, 

regardless of whether the medium is homogeneous or not. In an inhomogeneous 
medium 77 will be a function r. For the moment we shall assume that it is a 
constant. Solving for cj>' , we obtain the equation of the orbit as 

, , r a dr 



r \J {w f - r 2 a 
Integration of (JBJ gives the trajectory 



(6) 



0-0 o = cos-M^), (7) 

where <fio is a constant of integration. Squaring both sides of and using J5J 
give the conservation of energy 

f 2 + r 2 4> 2 + 2U(r) = -A, 

provided the radius of the caustic |Keller, Rubinow 1960| is given by 

r a = r 2 cf) = const. (8) 

This is Kepler's law of equal areas in equal times, or the conservation of an- 
gular momentum, where r 2 tj) is the angular momentum (relative to unit mass). 
Introducing © into J5J| results in 

rft = V(vr) 2 -rl (9) 



V 



\Jr 2 - r 2 a /r] 



The eikonal is the integral over r of the Legendre transform of the integrand 
of Fermat's principle 



= I - r a {<p - 0o ) = yj (i]r) 2 - r 2 a - r a cos 1 ( - 

\1 



r a I tan 2 (f)d(f> = r a (tan0 - <f) > 0. (10) 
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The second line of l(TU)) expresses the eikonal as the difference between Q and 
the arc length along the caustic |Keller, Rubinow 1960| . The third line of l|10[l 
shows that the product of r a and (|10|> is the area of a circular cap whose base 
is a circular arc of radius r a , and whose peak is a distance r from the center 
of the circle |Sholande"r"l 952 . In other words, the addition of a circular cap to 
a body increases its area by r a S(r), while its perimeter is increased by 2S(r), 
which is strictly increasing and strictly convex on r a < r < oo. 

The second line of H1(J|) is Debye's asymptotic expression for the phase of the 
Bessel function in the periodic region r > r a |Uebye 1 988 . In a homogeneous 
refractive body of a constant index of refraction, a light source is situated a 
distance r from the center of a circle of radius r a , which is the caustic. The 
eikonal l|10fl consists of two rays: a half-line \Jr 2 — r 2 a from the source to the 
point of tangency to the circumference of a circle of radius r = r a , and a 
ray along the arc length, r a cos~ 1 (r a /r). This quantity is subtracted from the 
straight line segment because the ray is taken from the caustic to the source. For 
the reverse path, the signs of the two terms are exchanged. Apart from a term 
— j7T, expression (|10fl . when multiplied by ft, represents the phase going away 
from the caustic, while the ray going toward the caustic has the signs reversed 
in addition to a phase factor of Hence, the phase changes abruptly by an 

amount — \-k upon passing through the caustic. A phase jump of this magnitude 
is observed for a focal line due to the convergence of rays of a cylindrical wave 
jSommerfeld 1964j . 

In the shadow region r < r a , where the rays do not penetrate, the eikonal 
ifTHjl becomes completely imaginary 

5l(r)=z{r Q cosh(^)- v ^| — 7^). (11) 

Since the 'shadow' intensities vanish rapidly as k, or the distance from the caus- 
tic increases, they are usually ignored [Poston and Stewart 1 978 . However, the 
matching conditions between the periodic and exponential regions are of fun- 
damental importance in quantum mechanics because they furnish the quantum 
conditions |Lanczos 1961| . Hence we must give credence to the type of motion 
that occurs in the shadow zone. Interestingly enough, contradictions will arise 
with the laws of physics that govern the illuminated, or periodic, zone. 

For instance, by reinstating the velocity of light c, we have the inequality 
r<t> > c in the shadow zone. This is in contradiction with special relativity. 
However, in the spectrum of a medium of anomalous dispersion there can exist 
a region near the absorption line where the group velocity can be greater than 
c |Brillouin I960] , In this region, the group velocity no longer represents the 
velocity of the signal. Anomalous dispersion that results from strong absorption 
destroys the characteristic wavelength of propagation so that light can only be 
defined statistically. For electrodynamic guided waves [Brillouin 19601 p 149], 
and in quantum mechanics, the phase velocity is greater than c. The concept of 
a phase velocity in quantum mechanics was discarded when Schrodinger [1982] 
proved that the group velocity of the wave represented the particle velocity of 
the electrons. But this meant that individual waves could not be used for signal 
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transmission. Even if they could be there would still be a detection problem 
since no optical effect could propagate with a velocity greater than c. 

Boundary conditions in general relativity are usually stated by requiring 
spacetime to be asymptotically fiat, as in the case of the Schwarzschild met- 
ric. But, in rotating systems a cutoff must be introduced for otherwise dis- 
tances r > c/4> would make the time component of the metric tensor negative 
|Landau and Lifshitz 1975| . Such a restriction would limit the phase of the 
Bcsscl function to the periodic region, but it would not be in the form of an 
asymptotic boundary condition. Larger distances where the angular velocity 
becomes greater than the velocity of light would bring us within the caustic 
region, and does not lead to the conclusion that such a system cannot be made 
up of real bodies |Landau and Lifshitz 1975| . 

The eikonal (|llf) can also be derived from Fermat's principle, which now 
reads 



/ = 



J rj\J r 2 <ft 2 — r 2 dt = J 7]\/ r 2 (f)' 2 — 1 dr = extremum. (12) 
Following the same procedure as before, we find the trajectory 

0-00 = cosh" 1 fcj , (13) 

for a constant index of refraction, where we can always arrange for O — by 
suitably fixing the initial point for the measurement of the arc length. The 
extremum (|12|) is just the distance / = — \Jr 2 — (r/r) 2 = —r a tanhc/f). 

Using the canonical parameterization for which r a — 1, the arc length s = 
sinh <j) enables the profile curve to be written as 

= {g(s),h(s)) = ( sinh" 1 s - / , J— 

where g(s) measures the distance along the axis of revolution, and h(s) measures 
the distance from the axis of revolution. The parameterization of the surface of 
revolution is 



x(s, 9) = ( sinh 



cos W sin ( 



\/l + s 2 Vl + s 2 Vl + s- 



where is the angle through which the profile curve has been rotated. The 
element of arc length on the surface of revolution is 

da 2 = dr 2 + r 2 sech 2 (tanh 2 6d6 2 + d9 2 ) . 



The eikonal Ijllfl is just the distance along the <fi axis, 



Sl(s) =sinlr 1 s- 



.s 



which in terms of 4>, 



S] 3 {(j))=r a { tanh 2 <p d<p = r a {(j>- tanh <f>). (14) 
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is a tractrix, having a constant negative curvature, K = — 1/r 2 . The tangent 
to the tractrix which intersects the x-axis always has the constant value r a . 
The distance from the origin to the point of tangency is r a (f>. The point on 
the tractrix which has a tangent intercepting the cc-axis is located a distance 
r a {<l> — tanh</>) along the x-axis. This is precisely the eikonal (|14f> . Hence, in 
the periodic domain the eikonal is half the increase in the perimeter due to 
the addition of a spherical cap to a body, while in the exponential domain, the 
eikonal is a tractrix, which is the involute of a catenary unwinding from its 
lowest point. 

The solutions to Helmholtz's equation are now exponentials 

^(0 = ^=6^(1 + 0(0)- 

The general solution has two free constants, B + and B~ , associated with the 
± signs in the exponent. The exponential growing solution is unphysical, and 
consequently, it must vanish. The matching condition between the constants A 
and 1} in (0J and B + is Lanczo s"l961| 

B + = \Acos{d- ±tt), 

and if it is to vanish, 

■& = -\lT±kTT, (15) 

.1 

where k is an arbitrary integer. This gives rise to the phase jump e*4 7r as the 
system passes through the caustic. 



3 The Mechanics of Diffraction Phenomena 

Formulas JHJl and JHJ are familiar from Kepler's theory, but depend on the form 
of the index of refraction. In a repulsive held U is positive, or in no held at all, 
the realness of the index of refraction J2J| requires A < so that the trajectory 
is open, extending between r = 00 and a minimum value of r, where the velocity 
vanishes. 

In an attractive held of force, A can be of either sign. The unperturbed 
motion will be the result of some playoff between the gravitational potential, 
U = lZ/2r and the centrifugal potential, r 2 /2r 2 . Introducing the former into 
the expression for the index of refraction © , the eikonal (|1U|I becomes 



y/-,4r 2 +Tlr-rl 
S{r,r a ) = / dr 



\J —Ar 2 + IZr — r 2 — r a cos 11 
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where the eccentricity 



4Ar 2 

(17) 



TV ' 

the semi-latus rectum 

q = 2 -£ = a(l-S), (18) 

and a = 1Z/2\A\, which is the semi-major axis if e < 1. The equation for the 
orbit is 

r = i " 1 ( 19 ) 

1 + e cos <j) ' ' 

which is an ellipse for e < 1 (A > 0), or a hyperbola for e > 1 (A < 0). In the 
periodic case, <f> is known as the 'true anomaly', and u is the 'eccentric anomaly'. 
The latter satisfies the equation 

r = a(l — ecosu) 

in the periodic case, A > 0, while it satisfies 

r = a (ecoshu — 1) , 

in the hyperbolic case A < 0. 
The velocity 



r = ±^- Al * +nr -% (20) 
r 

will vanish at 

r± = o(l±e). (21) 

r undergoes a libration, where the turning points are the aphelion (r + ) and the 
perihelion (r_). In the hyperbolic case, there will be only one point where the 
velocity H20JI vanishes, and that is at the closest distance of approach 

r min = a(e - 1). (22) 

A Hankel function of type j, with argument ner, and order n(q — r) has the 
integral representation 

^ fo _ r1 H = - [ e^r^-il-r)*) #j j = lj2 . (23) 
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The contours Cj are the paths from — 7r + ioo to — ioo and — ioo to tt + ioo. We 
seek an asymptotic expansion of 1)23(1 as k — > oo. On the basis of the definitions 
of the contours Cj , the only critical points of 

W(<t») = erBin<j>-(q-r)<l> (24) 

are simple saddle points. The necessary condition for simple saddle points, 
W((j>±) = 0, gives equation lfT^|l . where 4>+ = — (/>- w hh < cj) + < tt and — n < 
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4>- < 0. These saddle points will be simple because W"(<f>±) — —ersm(f)± ^ 0. 
Evaluating i|24[) at the saddle points leads to 

W(4>±) = ± y/-Ar 2 +Kr- r 2 a - (q - r) cos" 1 (^rj J ■ (25) 

As e — > co, implying that A < and the path is an hyperbola, (|25() becomes 
proportional to the free eikonal (|10fl . The derivative of 1251) with respect to r is 
proportional to the velocity (|2f)|> in this limit. We will have need of this limit in 
the next section when dealing with the bending of light by a massive body. 

If q — r > er, the order of the Hankel function 123fl is greater than its 
argument. The saddle points (f>± now lie on the imaginary axis 

4>± = iicosh -1 ( 

\ er 

and the phase becomes completely imaginary 

W\(t>±) = ±i | y/(q - r) 2 - e 2 r 2 - (q - r) cosh" 1 \^rj } ■ 

This corresponds to the shadow region where the rays do not penetrate. Since 
the contributions from (f> — are exponentially smaller than those from <j) = 
(/>+, the former can be neglected; hence, to leading order 

r(i) / 2 



exp |k Uq - r)cosh" 1 ((g/r - l)/e) - (q/r a )^Ar 2 ~TZr + r 2 a + {-l)H±i^ 



$jAr 2 -Kr + rl 

for j = 1,2. The change of signs under the square root sign imply that the grav- 
itational potential is now repulsive, while the centrifugal potential is attractive. 
In the periodic case, the shadow region lies beyond the limits of libration (|21|) , 
whereas in the aperiodic case, it is a circle whose radius is smaller than (|22() . 
Mechanically, these regions are classically inaccessible while optically, they can 
be penetrated, like quantum mechanical tunneling, and the breakdown of struc- 
tural stability in catastrophe theory. Near a caustic of the free phase Besscl 
function, it would be approximated by a cubic resulting in an Airy function. 
The catastrophe is known as a fold catastrophe Poston and Stewart 1978| . 

Finally, in the case A — 0, corresponding to a parabolic orbit (e = 1), the 
eikonal is 



S(r,r a ) = J V Ur ^ r * dr = 2y /TZr - r 2 a - 2r a cos" 1 



'Kr, 

2r a (tan0- <f>) > 0, (26) 



where 4> — cos 1 \ r a /yRr\ . The eikonal Ij26|l is similar to the unbounded 

motion, periodic motion of (|10f) . with the exception that the caustic radius r a 
is magnified by a factor r a /lZ. 
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4 Deflection of Light 



If U is the gravitational potential, the interaction is attractive, and the orbit 
curves towards the center of force. Gravity makes the medium optically more 
dense in the vicinity of the sun than far away from it. As a result, light waves 
will be bent around the sun rather than being straight. The effect was originally 
predicted by Soldncr in 1801, and rederived by Einstein in 1911 on the basis of 
the Doppler effect and Newtonian theory. Specifically, Einstein considered the 
slowing down of light in a gravitational field [Einstein 1952) . General relativity 
predicts a value twice as great, as Einstein showed in 1915 Einste in 1952) . It 
is therefore concluded that Newtonian theory is only an approximation, valid 
at speeds small compared to that of light, and the full theory must be used 
to calculate the actual deflection of a light beam traveling through the sun's 
gravitational field. The value obtained from Newtonian theory is only half that 
predicted from general relativity If instead the interaction were quadrupolc in 
nature, the phase of the Bessel function would belong to the periodic domain, 
and thus correspond to a diffraction phenomenon. 

The directions of the asymptotes of the hyperbola (|19fl are determined from 
the condition r = oo which gives 

A0 = cos -1 (-- J = 7r — cos -1 ( - j = \tx + skr 1 (- 



(27) 

Since only the rest energy is involved, A = —1, and e = 2r a /lZ for r a ^S> 1Z, the 
total deflection will be twice l|27l) or 



2A0 = 7r + 2sin~ 1 



K 

2^ 



A straight ray corresponds to 2A0 = tt. This is comparable to Coulomb scatter- 
ing, except that the trajectory is concave toward the origin where the massive 
body is located. For small angles of deflection the angle between the two asymp- 
totes differs from w by 6 = lZ/r a , which is the Newtonian result, but differs 
from the actual gravitational value by a factor of two. This is used to high- 
light the shortcomings of Newtonian theory and the need of general relativity 
|M0ller 1952] . 

M0ller [1952] splits the deflection of light into two effects: a velocity of 
light that varies as 1/ yl — TZ/r, and the non-Euclidean character of the spatial 
geometry. The first effect gives the radial velocity 



-f = H3>, (28) 



1 - TZ/r r 2 

to first order in TZ/r, while the second effect has a radial velocity 



1---^ + ^. (29) 
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Both give a contribution that is half the general relativistic value. Equation (|29|) . 
which arises from the Schwarzschild exterior metric tensor, has the wrong sign 
of gravitational field, ft appears repulsive like that of the centrifugal potential. 
Rather, we shall now show that the deflection of light by a massive body is 
accountable by the last term in 1|29|) . which is a quadrupole interaction, without 
the incorrect gravitational potential. 

When the light ray enters the scattering potential, it enters a region of 
changing index of refraction, where it is both refracted and diffracted. This is 
determined by the phase of the wave function. Since the gravitational potential 
is not responsible for the bending of light, we have to look to a higher order 
interaction. Expanding the gravitational energy of the central field of force in 
inverse powers of r we have 

2L/=-^(l + Cl ^ +C2 (^) 2 + ...), (30) 

where the Cj are coupling constants, and r a is the characteristic length. This 
is analogous to the Coulomb case, where r a would represent the radius of the 
first Bohr orbit [Sommerfeld 1950j . The next lowest interaction, ci, would be a 
mass dipole which, however, cannot oscillate so that its center of mass cannot 
accelerate. Hence, we set ci = 0, and consider the quadrupole term, Q/r 3 , 
where Q = Mr\ is the quadrupole moment of a system consisting of the sun 
and a photon in the direction perpendicular to the asymptote of the photon 
trajectory and the parallel line passing through the center of the sun. We shall 
consider the coupling constant c-i to be of order unity. 

It proves convenient |M0ller 19521 p 354] to introduce the new variable a = 
r a p\J\ — Hp into the equation of the trajectory, 

(31) 




Neglecting the small term Tip 3 , the equation of the orbit is Q, which is equiv- 
alent to 

j. dS 

4> + t; — = 4>o = const., 
or a 

where S is the eikonal of the Bessel function, (|lt)(l . For a constant index of 
refraction, the trajectory r — r a jr\ sin 4>, obtained by setting 4>q — \-k, is a 
straight line which passes through the origin at a distance r a /ri when cf> = 
and goes to infinity again for <p — > 7r. 

The exact equation (|31|l may be written as |M0ller 1952| p 354] 



dp 



-Vl-^ 2 , (32) 



where a = r a p^l — Hp. Since Hp is a small quantity, we obtain the approxi- 
mations 

a = r a p (1 - \Hp) r a p = a (1 + \Tlp) = a (l + -^-a 
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to first order, and consequently, r a dp = (1 + TZa/r a )da. Introducing these 
approximations into i|32|) . and integrating from to the distance of closest ap- 
proach, (To, we get 

f° da(l + Ka/r a ) K , s 

A0 = 7T + 6» = 2/ — v 1 al =7T + 2— . (33) 







In the derivation of H33J) we have used the fact that the distance of closest ap- 
proach is determined by the vanishing of the integrand in the denominator, and, 
hence, ctq = 1. Consequently, the deflection 9 = 21Z/r a is twice that obtained 
by treating the interaction through a Newtonian potential. The quadrupole 
interaction is introduced by the index of refraction, while the geometry is deter- 
mined by the flat metric, proportional to Q. Can we associate a wave function 
to such a process? 

For the quadrupole interaction we have a wave function whose eikonal and 
amplitude are approximately given by 

S Q (r) = ^-rl (l-^-r.cos" 1 fefi^jj , (34) 

and 

respectively, since IZ/r is a very small quantity. The square of the amplitude, 
(|35l) . is an intensity which measures the 'geometric divergence' of the wave field 
|Babic, Buldyrev 1991| . The larger it is, the more the rays in the neighbor- 
hood of some fixed ray will diverge from each other. In the absence of the 
quadrupole interaction the geometric divergence becomes infinite on the caustic 
r = r a , which is the case of the ordinary Bessel function. In the presence of the 
quadrupole term, the geometric divergence of the ray field becomes infinite at 
a distance smaller than r a . The equation of the trajectory, 



sin <p V r 

shows that the trajectory becomes infinite as <f> — > and ir, while the distance 
of closest approach r is less than the caustic radius r a . Consequently, as the 
result of the quadrupole interaction the distance from the scattering center 
to the caustic has been reduced. In the limit as r —> tq, the form of the wave 
function (@J is no longer valid, and a uniform asymptotic expansion of the Bessel 
function can be given in terms of Airy functions |Ludwig 1970| . 
The differential cross section of the quadrupole interaction is 



A(6>) = 2nr 



dr„ 



d0 



2 



which displays the characteristic scattering property of the cross section to ap- 
proach infinity as the angle of deflection approaches zero. However, since the 
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rest energy is involved, there will be no rapid decrease of 9 with the kinetic 
energy, as is usually the case in particle scattering. 



5 The Perihelion Shift 

Thus far we have treated open trajectories A < 0. In the case of Mercury, we 
are dealing with a closed orbit, and hence negative total energy A > 0, which 
is not exactly elliptical due to the advance of the perihelion. This precession 
amounts to 43" in a century, and was known to astronomers as early as 1860. 
However, its explanation had to await the general theory of relativity in 1915. 

Taking into account both the gravitational potential and the quadrupole 
interaction, the index of refraction is [vid. (|30|l with c\ = and C2 = 1] 




n -- 

In the unperturbed state, where the quadrupole interaction is absent, the gravi- 
tational potential must be large enough so that the eccentricity is real [vid. 117(1 
below] . The eikonal is 



Sp{r , r a ) = I vVEEll dr= [ V-Ari+Kr-rl+Krl/r ^ 



A closed trajectory will result from a dynamic balance between gravitational 
and centrifugal forces. 

In the presence of the quadrupole interaction, the angular momentum will 
exist only for radii r > ^1Z. Introducing this fact into through the trans- 
formation, r + §7?. = r' , and retaining terms that are at most quadratic in the 
Schwarzschild radius, give 

Sp(r,r a ) : 

where for brevity we have dropped the prime on r. Expanding the integrand in 
powers of the small correction terms results in 




S P (r,r a ) ~ 4 0) (r,r a ) + i^ 2 J 



dr 



^J-Ar 2 +TZr - r\ 



-Tlr't I dr = + ■■■ (37) 

V — Ar 2 + TZr — r 2 



where the unperturbed eikonal Sp\r,r a ) is given by H16J1 . The unperturbed 

:(0), 



orbit, —Sp/dr a — <p is given by the ellipse (|TU|l since e < 1. 



'p 



The change in Sp^ over one complete orbit is 



or a 



14 



As r goes through one libration, the true anomaly <j> increases by 2ir, and there 
is no perihelion advance. However, taking into consideration the first order 
correction term in (|37fl , which can be written as 

AS ( p>(r,r a ) = -f = |tt— , 

r a dr a r a 

and differentiating it with respect to r a , give the first order correction term 

dAS™ _ 3 TZ 2 
dr a 2 r 2 a 

Introducing the semi-latus rectum (|18|l results in 



= — «r- = ( 38 ) 



A0! = 3tt ^ (39) 
a(l — e z ) 

which is the general relativistic result for the perihelion shift. For Mercury, the 
rotation of the perihelion per revolution amounts to 0.104". The dimensionless 

3 

energy constant A = 2.59x 10 , and the mean motion oj = r a c/ab — 2cAz /TZ — 
8.34 x 10 -7 s _1 , where we have reinstated the speed of light c, and b = r a j \[A is 
the semi-minor axis. The period of the motion of Mercury is T = 2h/uj = 87.25 
days, which is close to the actual value of 88 days. The frequency of rotation of 
the perihelion will be Au = wA0i = 4.25 x 10~ 13 s^ 1 . 

There are five experimental confirmations of general relativity: time delay 
in radar sounding, deflection of light, perihelion advance, spectral shift, and the 
geodesic effect. We have shown that three of these effects can be treated as 
diffraction phenomena on the basis of Fermat's principle and the modification 
of the phase of a Bessel function in the short-wavelength limit. The spectral 
shift can also be derived from energy conservation by assuming that a photon's 
energy has both inertial and gravitational mass. This depends in an essential 
way on the equivalence principle, and not upon general relativity. 
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